Noname manuscript No. 

(will be inserted by the editor) 



Classical model of confinement 

Yu. P. Goncharov • N. E. Firsova 



Received: date / Accepted: date 



Abstract The confinement mechanism proposed earlier and then applied successfully 
to meson spectroscopy by one of the authors is interpreted in classical terms. For this 
aim the unique solution of the Maxwell equations, an analog of the corresponding 
unique solution of the SU(3)-Yang-Mills equations describing linear confinement in 
quantum chromodynamics, is used. Motion of a charged particle is studied in the field 
representing magnetic part of the mentioned solution and it is shown that one deals 
with the full classical confinement of the charged particle in such a field: under any 
initial conditions the particle motion is accomplished within a finite region of space so 
that the particle trajectory is near magnetic field lines while the latter are compact 
manifolds (circles). An asymptotical expansion for the trajectory form in the strong 
field limit is adduced. The possible application of the obtained results in thermonuclear 
plasma physics is also shortly outlined. 

Keywords Quantum chromodynamics • Confinement • Thermonuclear plasma 
physics 



1 Introduction 



In Refs. 1,2,3] for the Dirac- Yang-Mills system derived from QCD-Lagrangian an 
unique family of compatible nonperturbative solutions was found and explored, which 
could pretend to decsribing confinement of two quarks. The successful applications 
of the family to the description of both the heavy quarkonia spectra [4jj5j[6l[7j and 
a number of properties of pions, kaons, r\ and r/' -meson 8 9, 10 11,12 showed that 
the confinement mechanism is qualitatively the same for both light mesons and heavy 
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quarkonia and it is mainly governed by the magnetic colour field linear in r (distance 
between quarks) which represents a part of the mentioned unique family of solutions 
and, in its turn, the part is a solution of the SU(3)-Yang-Mills equations. 

As has been emphasized, however, as far back as in Refs. 2,3 , the similar unique 
confining solutions exist for any compact semisimple and non-semisimple Lie groups, 
in particular, for SU(iV)-groups with N > 2 and for U(iV)-groups with N > 1, i.e. it 
holds true also for classical electrodynamics with group U(l) and Maxwell equations. 
Under this situation, as was pointed out in Refs. 2 ,3 , there is an interesting possibility 
of indirect experimental verification of the confinement mechanism under discussion. 
Indeed the confining solutions of Maxwell equations for classical electrodynamics point 
out the confinement phase could be in electrodynamics as well. Though there exist 
no elementary charged particles generating a constant magnetic field linear in r, the 
distance from particle, after all, if it could generate this elecromagnetic field configura- 
tion in laboratory then one might study motion of the charged particles in that field. 
The confining properties of the mentioned field should be displayed at classical level 
too but the exact behaviour of particles in this field requires certain analysis of the 
corresponding classical equations of motion. 

The aim of the present paper is to some degree to realize the above program on 
studying motion of the charged particles in the mentioned confining electromagnetic 
field. 

Section 2 contains preliminaries necessary to pose the problem: information on the 
confining solutions of the Yang-Mills and Maxwell equations and on the miscellaneous 
forms of the motion equations for a charged particle in the confining magnetic field 
when considering it with using different curvilinear coordinates. Section 3 is devoted 
to the general conclusions of a qualitative character concerning behaviour of a charged 
particle in the magnetic field under discussion. In the strong field limit Section 4 gives 
asymptotical expansions for the spherical coordinates of a particle when its moving in 
the field under consideration while Section 5 contains numerical estimates and Section 
6 is devoted to discussion and concluding remarks. 

Appendix A is devoted to the formulation of vector analysis on a region J? in R 3 
which is most convenient, especially while working with using the arbitrary curvilinear 
coordinates so the mentioned formulation is employed throughout the paper. At last, 
Appendix B supplements Section 2 with a proof of the uniqueness theorem from that 
Section in the case of U(l)-group (Maxwell equations). 

Also throughout the paper we employ the Heaviside-Lorentz system of units with 
h — c = 1 and also with the Boltzmann constant k = 1, unless explicitly stated 
otherwise. When calculating we apply the relations 1 GeV^ 1 ~ 0.1973269679 fm , 

1 s _1 « 0.658211915 x 10~ 24 GeV , 1 V/m « 0.2309956375 x 10~ 23 GeV 2 , 1 T = 
4tt x 10~ 7 H/m x 1 A/m « 0.6925075988 x 10~ 15 GeV 2 . 

2 Preliminaries 

2.1 The confining solutions of SU(3)-Yang-Mills and Maxwell equations 

As was mentioned above, our study is motivated by the confinement mechanism pro- 
posed earlier by one of the authors and based on the unique family of compatible non- 
perturbative solutions for the Dirac- Yang-Mills system (derived from QCD-Lagrangian) 
studied at the whole length in Refs. [1,2,3 . 
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One part of the mentioned family is presented by the unique nonperturbative confin- 
ing solution of the SU(3)-Yang-Mills equations for gluonic field A = A^dx^ = A^Xadx^ 
(A are the known Gell-Mann matrices, fi = t, r, #, <p, a = 1, 8). To specify the ques- 
tion, let us note that in general the Yang-Mills equations on a manifold M can be 
written as 

d* F = g(*F A A - A A *F) , (1) 

where the curvature matrix (field strentgh) F = dA + gA AA = F^ v X a dx^ A dx v with 
exterior differential d and the Cartan's (exterior) product A, while * means the Hodge 
star operator conforming to a metric on manifold under consideration, g is a gauge 
coupling constant. 

The most important case of M is Minkowski spacetime and we are interested in the 
confining solutions A of the SU(3)-Yang-Mills equations. The confining solutions were 
defined in Ref. [I] as the spherically symmetric solutions of the Yang-Mills equations 
(1) containing only the components of the SU(3)-field which are Coulomb-like or linear 
in r. Additionally we impose the Lorentz condition on the sought solutions. The latter 
condition is necessary for quantizing the gauge fields consistently within the framework 
of perturbation theory (see, e. g. Ref. [13]), so we should impose the given condition 
that can be written in the form div(yl) = 0, where the divergence of the Lie algebra 
valued 1-form A = A^dx 1 * = A^Xadx^ is defined by the relation (see, e. g., Refs. [191 

W) 

div(A) = ^du(VSg^A v ) . (2) 

It should be emphasized that, from the physical point of view, the Lorentz condition 
reflects the fact of transversality for gluons that arise as quanta of SU(3)-Yang-Mills 
field when quantizing the latter (see, e. g., Ref. |13j). 

Under the circumstances, the unique nonperturbative confining solution of the 
SU(3)-Yang-Mills equations looks as follows 

A 3 4- 1 A s - ai 4- A A 3 4- 1 A 8 - 0,2 4- A 

A t+ —A t ---+A 1 ,-A t+ —A t ---+A 2 , 
9 4 = £i±£2_ (Ai+A2)) 



V3 

A% + —A% = bir + Bx , -A% + -^A% = b 2 r + B 2 
v3 v 3 



--j=A% = -{b 1 +b 2 )r-{B 1 +B 2 ) (3) 

with the real constants aj , Aj , bj , Bj parametrizing the family. As has been repeatedly 
discussed by us earlier (see, e. g., Refs. [2,3. and below), from the above form it is clear 
that the solution (3) is a configuration describing the electric Coulomb-like colour field 
(components A^ ,S ) and the magnetic colour field linear in r (components ^4^' 8 ) and we 
wrote down the solution (3) in the combinations that are just needed further to insert 
into the corresponding Dirac equation (for more details see Refs. [T|[2j[3] ) . 

The word unique should be understood in the strict mathematical sense. In fact in 
Ref. [2] the following theorem was proved (see also Appendix B): 

The unique exact spherically symmetric (nonperturbative) confining solutions (de- 
pending only on r and r^ 1 ) of SU(3)-Yang-Mills equations in Minkowski spacetime 
consist of the family of (3). 
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It should be noted that solution (3) was found early in Ref. [T] but its uniqueness 
was proved just in Ref. [2] (see also Ref. [3]). Besides, in Ref. [5] it was shown that the 
above unique confining solutions (3) satisfy the so-called Wilson confinement criterion 
14, 15]. Up to now nobody contested this result so if we want to describe interaction 
between quarks by spherically symmetric confining SU(3)-fields then they can be only 
those from the above theorem. 

Now one should say that the similar unique confining solutions exist for all semisim- 
ple and non-semisimple compact Lie groups, in particular, for SU(JV) with N > 2 and 
U(7V) with N > 1 0|3]. Explicit form of solutions, e.g., for SU(iV) with N = 2,4 can 
be found in Ref. [3] but it should be emphasized that components linear in r always 
represent the magnetic (colour) field in all the mentioned solutions. Within the present 
paper we are especially interested in the U(l)-case (electrodynamics) and a proof of the 
above uniqueness theorem for that situation is adduced in Appendix B for inquiring. 

Under this situation the Yang-Mills equations (1) turn into the second pair of 
Maxwell equations 

d * F = (4) 

with F — dA, A = A^dx^ . As is discussed in Appendix B, in the spherically symmetric 
case the equations (4) are equivalent to 

d r (r 2 d r A f ) =0, drA v = 0, (5) 

with At = At(r), Aip — Atp(r) and we write down the unique solutions of (5) as 



A t = -+A,A v = br + B (6) 



r 

with some constants a, b, A, B parametrizing solutions. 

To interpret solutions (6) in the more habitual physical terms let us pass on to 
Cartesian coordinates employing the relations 

<p — arctan(j//s), dtp = i^~dx + ~^~dy (7) 

which entails 

. , /, „s , (br + B)y , (br + B)x , , , 

A = Amdip — (br + B)dtp — — = tfdx + ^ '^-dy (8 

x- 2 +y £ x £ + y A 

and we conclude that the solutions (6) describe the combination of the electric Coulomb 
field with potential $ = At and the constant magnetic field with the vector-potential 
(8) which can be written as (using isomorphism dx i, dy j, dz <=> k, see 
Appendix A) 

A = ^i + ^j + ^k = - (6 ^i + (6 ^ j , (9) 

x + y x + y 

which is linear in r in spherical coordinates. Let us compute 3-dimensional divergence 
divA with the help of 3-dimensional Hodge star operator (see Appendix A). In spherical 
coordinates we have [see (A. 12) and (A. 7)] 

divA = *(d * A) = *d * [(br + B)dip] = *d ( bT + f dr A d# ) = . 
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Then eqs. (4) in Cartesian coordinates take the form 

A$ = 0, rot rot A = grad divA - AA = -AA = (10) 

with the Laplace operator A = d 2 + dy + d 2 . Also it is easy to check that the solution 
under consideration satisfies the 4-dimensional Lorentz condition &\v(A) = 0, where 
4-dimensional divergence is defined by (2). 

Finally, as is shown in Appendix A, the corresponding strength of magnetic field is 

, , . •. b b f xz yz , \ 

H = rot A = *(dA) = do = ^dx H ^f- ^dy — dz 

sinw r \x* + y 4 x z + y z J 

b I xz yz . , \ , . 

2^ 2 1 + 2T 2 J~ k ' 11 

r \x z + y z x z + y z J 
respectively, in spherical and Cartesian coordinates. 



2.2 Singularities of solutions 

As is seen from (8) and (11), the magnetic field under exploration has the singularities 
on the z-axis so its formal mathematical definition domain is the manifold R 3 \{z} 
with the z-axis discarded rather than the manifold M 3 . Singularities of such a kind 
are of mathematical nature and appear when trying to write a concrete macroscopic 
physical field in an analytical form. Physical origin of the given singularities is that 
some sources generating the field should be present on 2-axis. Another matter is that 
the field under consideration may probably be modelled by miscellaneous ways. For 
the sake of completeness, one of possible physical realization will be considered in Sec. 
6. In theoretical considerations within classical approach one should segregate from a 
concrete realization of one or another macroscopic electromagnetic field and consider 
them to be given on their natural mathematical definition domains. 

At the quantum level, however, treatment of singularities may be different from 
classical one. In particular, in the case of gluonic field (3) the problem of singularity 
along z-axis of magnetic part for solution (3) can be resolved by that quarks may emit 
gluons outside of some cone $ = t?o so singularity along 2-axis plays no role (for more 
details see Ref. 112; and estimates for $rj in pions and kaons therein). 



2.3 Equations of motion for a charged particle in the confining magnetic field 

As was mentioned in Section 1, at quantum level the confinement of quarks is ba- 
sically governed by the magnetic (colour) part (linear in r) of solution (3), as has 
been discussed in Refs. [4,5,6,7,8,9,10,11,12;. In the present paper we would like, at 
classical level, to explore the behaviour of a charged particle moving in the confining 
magnetic field (11). Accordingly, we need to study classical equations of motion for 
such a particle. As is known (see, e.g., Ref. |16p. those equations are obtained from 
Lagrangian 

L = -rn\J 1 - v 2 + qAv, (12) 

where q and m are, respectively, charge and mass of a particle while the form of both 
the velocity square v 2 = g^v^Vv and the scalar product Av = g^ v Ap,v v depends 
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on choice of curvilinear coordinates. Then the sought equations are derived from (12) 
according to the standard prescription of Lagrange approach as 



£ ( dL \ - ( dL 

dtyd&J \ d Qi 



0,i = 1,2,3, (13) 



where Qi are the chosen coordinates and the dot signifies differentiation with respect 
to t. For our purposes the equations of motion will be useful in both spherical and 
Cartesian coordinates. One can note that v 2 is conserved |16j when moving in a constant 
magnetic field, i.e., the full energy E — m/yl — v 2 of (relativistic) particle is also 
conserved. In the case of spherical coordinates we have v = r dr+r 2 ft dft+r 2 sin 2 ftip dip, 
v 2 = g^v^Vu = r 2 + r 2 ft 2 + r 2 sin 2 flip 2 = v 2 = const, Av = g w ' A v v v = (br + B)ip 
with A from (8) and in accordance with (13) we obtain 

a(r — r sin 2 flip 2 — rft 2 ) = <p , (14) 



£ 

dt 



^r 2 ft^j — r 2 (p 2 sin ft cos ft = , (15) 



jU— yr ip sin ftj — —f (16) 
with dimensionless parameter a = Ej (qb). In the case of Cartesian coordinates we have 

• j , • j i ■ j 2 -2, -2, -2 2 , A • (br+B)y , ■ (br+B)x 

v = xdx + ydy + zdz, v =x + y +z = v Q = const, Av = -x \. 2+y .i y + y \- 2+y i 
with A from (8)-(9) and (13) gives rise to 

fj,x=l(y + z y * ) , (17) 
r \ x z + y z J 

ny = --(x + z * Z 2 ) , (18) 
r \ x z + y z J 

nz = — 2Y ~ 2/*) • ( 19 ) 
r(a; z + y z ) 

Also we should add the initial conditions to (14)-(16) and (17)-(19). Namely, putting 
an initial moment of time to — for simplicity, we have, respectively, 

r(0) = r , ft(0) = ft Q , <p(0) = po, f (0) = r ,ft(0) = tf , 0(0) = pq . (20) 

or 

x(0) = x , J/(0) = 2/o, z(0) = 20, i(0) = x , y(0) = y , i(0) = i . (21) 



3 General considerations 

3.1 Magnetic field lines 

Let us above all find out how the magnetic field lines look for the field of (11). According 
to a general prescription (see, e.g., Ref. [21]) we should determine integral curves for 
differential system 

dx dy dz 



H x H y H z 



(22) 
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Fig. 1 Magnetic field lines of the confining magnetic field 



which can be made if finding the first integrals for it, i.e., such functions tp that satisfy 
the partial differential equation 

Then, as is not complicated to check, the system (22) has two independent first inte- 
grals, namely y = C\x, x 2 +y 2 + z 2 = C 2 with constants C\$,- That is, integral surfaces 
are planes and spheres and, as a result, integral curves are circles. For example, the 
plane y = is integral surface and equations of field lines are x 2 + z 2 = C 2 (see Fig. 
1). 



3.2 The confining properties 

We can note that 

d 2 d d 

-^r = — (2rr) = 2 — (xx + yy + zz) = 2[v + (xx + yy + zz)] (24) 

with v 2 — x 2 +y 2 + z 2 . Multiplying (17), (18), (19) by x, y, z, respectively, and adding 
the results, we get 

l i(xx + yy + zz)= xi ~ yX (l + z A , n = - . (25) 
r \ x z +y z J qb 

To calculate xy — yx we notice that ^ [x y — yx) = xy — y'x and replacing x, y according 
to (17) and (18), conformably, we shall have [lixy — yx) = — \.{xx + yy + zz) = —r 
wherefrom 

H(xy - yx) = -r + A , (26) 

where a constant Aq can be found from initial conditions (21) when considering (26) 
at t = so A = n(x yo - yo^o) + ro and 

(i(xy - yx) = - \/x 2 + y 2 + z 2 + fi(x Q y - y Q x ) + r . (27) 

We can consider Ao > which always holds true for the strong enough field when 
|6| — ¥ oo and, consequently, /x — > 0. Then, considering (27) on z-axis where x — y — 0, 
we obtain \z\ = Ao which signifies that particle trajectory can reach z-axis only at 
z = ±Aq. It should be recalled that according to Sec. 2 the z-axis is forbidden for 
motion of a particle so if for the particle \z\ — Ao at some moment then after it one 
should consider the motion to be finished and the particle vanished. Now, using (26), 
we can rewrite (25) as 

^(xx + yy + zz) = ( 1 + Z \ = ^ 4( ^ (1 + cot 2 tf) (28) 

r \ x z +y^J r 

with spherical coordinate #. At last, with the help of (28) we derive from (24) 

v 2 + ^_I(i + cot 2 tf) . (29) 



d 2 d 
d? r =dt {2rr) 
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Fig. 2 Behaviour of a charged particle in the confining magnetic field 



At r < Aq from here it follows 



v 2 + 



A 



" 7 -(l + cot 2 tf) 



dt > 



dt>0, 



(30) 



which signifies that r is increasing. But if r > 2Ao, i.e., r — > Aq, then Aq/v < 1/2 
and from (30) we gain 

r- A 



1 

r — - 
r 



-~2 + 
A 4 



Aq^ 
r/j 2 



provided that 



v < 



cor 



dt < 



1 

2Aq 



1 



< , (31) 



(32) 



y/2\fi\ s/2E ' 

i.e., r is decreasing. It should be emphasized that for sufficiently strong field (\b\ — > oo) 
the condition (32) will always be fulfilled for the given E. Besides, under this situation, 
Aq = /i(x'oyo — Vqxq) + r a ~ an d we can see that spherical coordinate r never tends 
to infinity and oscillates near the initial value ro. Inasmuch as, as said above, r = tq, 
f — ipQ is a magnetic field line, then we can say that the particle trajectory oscillates 
near the magnetic field line defined by initial conditions. In other words, we get the 
full confinement of charged particle in the magnetic field under discussion which is 
sketched out in Fig. 2. 



4 Asymptotical expansions 

To illustrate the general properties decribed in Section 3 it should be noted that the 
system (14)-(16) seems to be insoluble in an explicit form. But let us try to obtain an 
asymptotical solution of it in the form of expansions in the dimensionless parameter 
/i = E/(qb) in the strong field limit when b — > oo, i.e., /i — > 0. For this aim we can 
notice that the angle for Lagrangian L of (12) is the so-called cyclic coordinate, i.e. 
L does not depend on ip. Then in accordance with the Lagrange approach we have an 
integral of motion [see (13)] in the form 

9 L 2 2 2 2 

— = Er (p sin $ + q(br + B) = a v = Er a (pQ sin $o + q{brQ + B) — const , (33) 
oip 

with using the initial data of (20). From here it follows 

<i = \ . \ - , v = ro + M?Wo sin 2 # , (34) 
fir z sin v 

and it is not complicated to rewrite the system (14)-(16) in the form 



dt r A r 3 sin 2 ■& 

2 ds {y - r) 2 cost? 
dt r 1 sin 3 ■& 

r=p, (37) 

<>=^ (38) 

with s = r 2 i?. 



dp 2 s v(y — r) 
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4.1 Expressions for r and § 

Further we seek for r, i9, s, p in the form 

r = fo + cn+|i 2 f 2 + 0{i?) , = do + (idi + 0(a 2 ) , 

s = so + A* si + 0(/i 2 ) ,p = po + fJ.pi + 0(u 2 ) , (39) 

where f, fi, f2, #o, $1, so, si, po, Pi are some functions of t. Now, expanding the 
right-hand side of (35) in fx, we obtain 

M 2 (i> + i*Pi) + o(m 4 ) = r ° 3 (r ° i y + o(m) , 

rg sm z i?q 

so we should have fn = ro. 

Then the following terms of expansion for the right-hand side (35) give rise to 
relation 

2,_ . - \ , /->/■ 4n rgVo sin 2 $ - n , 2 . 

A* (Po + A* Pi) + 0(A» ) = _ 2 . 2 □ A* + 0(M ) , 

rg sin Vq 

which yields f\ = r 2 <^o sin 2 t?o = const. In this situation the new terms of expansion 
for the right-hand side (35) lead to 

a 2 (p Q + Hpi) + O^) = - 2 - 2 2 % ) » + 

and it should be 

Po = 4 ~ • (40) 

rg rg sin' i? 

Let us now pass on to the equation (36) where the conforming expansion with the help 
of (39) yields 

2/- , — \ I s—\ / 4x f|cosi?o 4 . 5n 
At (s + /iSi) + O(At )= 2 . 33 At +0(m)- 
rg sin° i? 

From here we have so = and, consequently, so = Co = const. Accordingly the 
equation (37) gives rise to fo + /ifi+0(/i 2 ) = Po + AtPl + 0(ai 2 ) which entails Pq = fo = 
0, pi = f 1 = since we have above obtained that fo = ro, fi = rgt^o sin 2 $0 = const. 
In the circumstances the relation (40) gives 



C 2 

f 2 = ^ sin 2 #0 • (41) 



ro 

At last, in a similar way from (38) we can obtain the relation 

^ + a^ 1+ O(p 2 ) = %+0(n), 
r o 

which entails #0 = % = and, as a result, $0 = + Cl with some constant C\. 
Then, taking into account (39) and (41), we finally have 

r = r + Atr^osin 2 ^ +At 2 — sin 2 (% + Cj VoO* 3 ),^ = % + C*i + O0*) . (42) 
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4.2 Expression for ip 

When searching for ip in the form <p = po + fi ip\ + 0(/j, 2 ) we shall, according to (34), 
obtain 

M£o + M<£i) + 0(a* 3 ) = — r2 2 a* 2 + 0(^ 3 ) , 

sin' $o 

■- - C 2 

wherefrom, with the help of (41), we get ip = 0, </5 1 = f- and, consequently, ipo = C2, 

r o 

C 2 

ipi = g-A + C3 with some constants C2 3. So finally 

r o 

= C 2 + (c 3 - ^ A* + 0(a* 2 ) • (43) 



4.3 Determination of constants 

At ji — > with taking (20), (42) and (43) into account we find Co ~ r #0i Ci ~ #0i 
C2 ~ ¥>o- This eventually leads to the final expressions 

2 2 2 3 " 2 2 3 

r « r + /ir <po sin i? + A* r #0 sin (i?o* + $0) + ), J? ~ #0 + #0 t + O(fx) , 

^~ p {) + (C 3 -r dlt)^ + O{pi 2 ), (44) 
that confirm the general considerations of Section 3 (see also Fig. 2) . 



5 Numerical estimates 

5.1 General estimates 

As is clear from (44), if we want a charged particle to be near magnetic field line 
r = ro, (p — <po defined by intitial conditions then we should impose the condition 
|/iro 2 y>o sin 2 #ol << r o which entails 

Atrwosin 2 ^! << 1 , A* = ^ • (45) 

If defining the physical components of velocity in spherical coordinates by equality 
« 2 = (v? h f + {vff + {vff = <T = r 2 +r 2 tf 2 +r 2 sin 2 V = «§ = r 2 +r 2 tf 2 + 
r\ sin 2 i?o^o ~ consi then we get v% — r, = ri9, = r sin i)(p and the condition 
(45) signifies that (vo)^ sin$o| < Im^oI << 1- 

But, obviously, _E = m/ \Jl — Vq so the condition \(xvo\ « 1 can be rewritten as 

« \b\ (46) 



v /(l-«g)47rJV 2 a em 



with g = JVe and the electromagnetic coupling constant a em = e 2 /(47r) « 1/137.036 
in the chosen system of units. 
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5.2 Deuteron in thermonuclear plasma 

For this case typical values ro ~ 1 m (see e.g. Ref. |17| ) and at temperature of plasma 
T ~ 0.8 625 x 10~ 2 MeV (10 s K) we have a mean thermal deuteron velocity vq ~ 
i/3r/m « 0.372 x 10~ 2 with the deuteron rest energy m = m p + m„ - 2.225 MeV « 
(938 + 939-2.225) MeV = 1874.775 MeV, iV = 1 and (46) yields \b\ » 23.0 MeV. Let 
us take |6| = 1 GeV and, for the sake of simplicity, put i?q = n/2. Then in accordance 
with (A. 17) the module of magnetic field strength near the particle trajectory will be 

H = — +L- ~ 0.197 x 10~ 15 GeV 2 ~ 0.351 T , (47) 
ro sin wo 

i.e., it is a quite accessible value under the laboratory conditions. It should be, however, 
noted that for the time t of order 1 s necessary to confine plasma before thermonuclear 
reaction starts [17] the angle •& can get increase 'dot = (vq)^ t/ro ~ vot/ro ~ 10 6 
according to (44), i.e. promptness of particle along the magnetic field line r = ro, 
tp = <po (see Fig. 2) will approximately be equal to Nq = 10 6 /(27r) ~ 1.77 x 10 5 , i.e., 
the particle will repeatedly cross z-axis which is impossible since the z-axis is forbidden 
for motion, as was mentioned in Sec. 2 and 3. But if (uo)$' 1 — > then Nq — > as well. 
We can draw the conclusion that deuteron rushing in to the field with (vq) 1 ^ 1 ~ will 
remain near its initial position during the time t = 1 s. 



5.3 Quarks in pions 

We may with minor reservations try applying the results obtained also to quarks within 
hadrons, e.g., within charged pions . In this case quarks are moving in the classical 
confining SU(3)-gluonic field (3) but they are quntum objects described by the wave 
functions - the modulo square integrable solutions of the Dirac equation in the field 
(3) (for more details see Refs. [8l l9lll2p . Let us, however, look at what the classical 
estimate (46) can give for quarks where, obviously, electric charge should be replaced 
by colour one and a e m by a s , the strong coupling constant. We can use the fact 8,9 
I12| that the colour magnetic field between quarks can be characterized by an effective 
colour strength H = 6/(rsini?) with b = \Jb\ + &1&2 + k 2 an d 61 2 from the solution 
(3) while r stands for the distance between quarks. Then (46) allow us to introduce 
quantity 

60 = . . (48) 

and for M-quark in 7r ± -mesons with m = m u ~ 2.25 MeV, vq ~ 0.99, a s ~ 0.485 
|12j we obtain 60 ~ 22.68 MeV which at the scale of pion ro ~ 0.672 fm entails 
H = 6 /(r sini7) ~ 0.666 x 10" 2 MeV 2 ~ 0.119 x 10 14 T when <d = tt/2. In reality, 
from quantum considerations [12; for vr^-mesons it follows b\ — 0.178915 GeV, 62 — 
-0.119290 GeV so b = ^Jb\ + b x b 2 + b\ w 0.157 GeV > b Q and at the scales of the 
meson under consideration we have [15] H ~ (10 15 - 10 16 ) T. As a result, classical 
estimate corresponds to the quantum considerations. 



6 Discussion and concluding remarks 



It is useful to compare our results with the well-known case of motion of a charged 
particle in the homogeneous magnetic field (see, e.g., Ref. [16p which is sketched out 
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Fig. 3 Motion of a charged particle in the homogeneous magnetic field 
Fig. 4 One possible physical realization of the confining magnetic field 

in Fig. 3. In the latter case the particle moves along helical curve with lead of helix h — 
2nmv cos a / {qH\/ 1 — v 2 ) and radius R = mv sin a /{qHVl-v 2 ). As a consequence, 
the homogeneous magnetic field does not give rise to the full confinement of the particle 
since the latter may go to infinity along the helical curve. Another matter is the case 
of the magnetic field (11). As we have seen above it provides the full confinement of 
any charged particle in case the field is strong enough: under any initial conditions 
the particle motion is accomplished within a finite region of space so that the particle 
trajectory is near magnetic field lines while the latter are compact manifolds (circles). 
This was explicitly demonstrated in Section 4 by obtaining asymptotical form of the 
motion under discussion. 

Taking into account such remarkable properties of the magnetic field in question we 
may hope that it should find a number of applications, in particular, in thermonuclear 
plasma physics where the problem of confinement of plasma during a sufficiently long 
time has so far not solved in a satisfactory way [17]. But for it one should explore the 
possible ways of modeling the field (11) in laboratory conditions which is seemingly 
not such a simple task. One of possible physical realization is sketched out in Fig. 4 
and is accomplished between two cone ferromagnetic pole pieces where the fact is used 
that magnetic field lines of ferromagnet are perpendicular to its surface. Of course, 
field lines inside the pole pieces (shown by dash lines) break off but if diameter D — > 
while the pole pieces are approaching, the whole construction tends to the the formal 
mathematical definition domain R 3 \{z}. In this realization it is clear why particles 
will leave the field: they will just be absorbed by pole pieces when moving along the 
field line. So one needs to fit parameters (in particular, the values of module H) of the 
whole construction in such a way that a particle could remain on trajectory for a long 
enough time t (e.g., for deuteron in thermonuclear plasma t~ 1 s according to Sec. 5) 
not reaching the pole pieces. 

Finally, as has been said in Section 1, the main motivation of writing the given 
paper was to interpret the mechanism of quark confinement proposed in Refs. 1,2,3 
in classical terms. As seems to us, we could to a certain degree do it. 

Acknowledgements Yuri Goncharov is grateful to Prof. F. Nasredinov from Experimental 
Physics Department of the Sankt-Petersburg State Polytechnical University for discussion 
regarding the confining properties of the magnetic field of (11) and for pointing out the possible 
physical realization of it. 



Appendix A 

In main body of the paper we employ the formulation of vector analysis on a region 
Q in R 3 from Ref. j!8| . To our mind such a formulation is most convenient especially 
while working with using the arbitrary curvilinear coordinates. The essence of that 
formulation is in systematical use of both the Hodge star operator * and the exterior 
differential d. 
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Hodge star operator on R and Minkowski spacetime 

Let M is a smooth manifold of dimension n so we denote an algebra of smooth functions 
on M as F(M). In a standard way the spaces as of smooth differential p- forms A P (M) 
(0 < p < n) are defined over M as modules over F(M). If a (pseudo)riemannian 
metric G = ds 2 = g^dx 11 (g) dx v is given on M in local coordinates x — (x^) then G 
can naturally be continued on spaces A P (M) by relation 

G(a,f3) = det{G(a i ,P j )} (A.l) 

for a = aiAtt2-Aa p , /? = fi\ A/?2---A/? p , where for 1-forms = of^ dx^ , /3j — fi^ dx v 
we have G(aj,/3j) = g^ v oijj} with the Cartan's wedge (exterior) product A. Under 
the circumstances the Hodge star operator *: A P (M) — > A n ~ p (M) is defined for any 
a € A p (M) by 

a /\(*a) = G(a,a)u) g (A. 2) 

with the volume n-form u g = a/| dct(g ^[dx 1 A ...dx n . For example, for 2-forms F = 
F^dx^ A dir^ we have 

F A *F = {g tia g vP - g^ g va )F^F a[j *f~5 dx 1 A dx 2 ■ ■ ■ A dx n , fj, < v, a < j3 (A.3) 

with 5 — | det(g M! ,)|. If s is the number of (-1) in a canonical presentation of quadratic 
form G then two the most important properties of * are 

* 2 = (-i)P("-f)+ s , (A.4) 

+ /2Q2) = /i(*«i) + /2(*a2) (A. 5) 

for any /i,/2 G F(M), a\,ai € yl p (M), i. e., * is a F(M)-linear operator. By virtue 
of (A. 5) for description of *-action in local coordinates it is enough to specify *-action 
on the basis elements of A P (M), i. e. on the forms dx 11 A dx 12 A ... A dx lp with i\ < 
12 < ••■ < ip whose number is equal to C p = ( w _!^)| p ; ■ 

The most important case of M in the given paper is a region Q in the Euclidean 
space R 3 with local Cartesian (x, y, z) or spherical (r, ■&, ip) coordinates. The metric is 
given by either ds = dx +dy + dz or ds = dr +r [dv + sin ddip ) , and we shall 
obtain the *-action on the basis differential forms according to (A. 2) in both the cases 
as 

*dx = dy A dz, *dy — —dx A dz, *dz = dx A dy, 

*(dx A dy) = dz, *(dx A dz) = — dy, *(dy A dz) = da;, *(d:r A dy A dz) = 1, (A. 6) 

2 1 

*dr = r sini9d# A dy>, *d$ = — sin $dr A dip, *dtp — -dr A d$, 

sm v 

*(dr A d-d) = sin-ddip, *(dr A dip) = ^d$, *(d$ A dtp) = - — -dr, 

sm » r z sm v 

*(dr A d-d A dip) = — -, (A. 7) 

H sin w 

so that on any p-form * 2 = 1, as should be in accordance with (A.4). 

Let us also adduce for inquiring the conforming relations for the case of cylindrical 

9 9 9 9 9 

coordinates p,ip,z, where x = p cos ip, y = p sin <p and metric ds = dp + p dip + dz . 
Then 

*dp = pdip A dz, *dip = dp Adz, *dz — pdp A dip, 

P 
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*{dp A dp) = -dz, *(dp A dz) — —pdp , *(dp A dz) = -dp, 
P ' P 

*(dp Adp A dz) = -, (A. 8) 

P 

where we can again see that * 2 = 1 on any p-form according to (A. 4). 

Also the important case of M is Minkowski spacetime with coordinates t, r, ft, p 
and metric ds = dt — dr — r [d-d + sin ft dtp ), where we have 

2 2 

*dt = r sin ftdr A dft A d</s, *dr = r sini9<f£ A dft A dip, 

*dft — — sini9df A dr A dip, *dp = - dt A dr A dft, 

sin i? 

*(dt A dr) = — r 2 sini?d$ A dip , *{dt A dft) = sin ftdr A dip , 

*(dt A dip) = —dr A dft , *(dr A dft) — sin ftdt A dip , 

sin ft 

*(dr A dip) = ^dt A dft , *(dft A dp) — —— - — -dt A dr , 

smft H sin ft 

*(dt A dr A dft) = sinftdip , *(dt A dr A dp) = -^dft, 

sine 

*(dt A dft A dp) = — rdr, *(dr A dft A dp) = — -dt, (A. 9) 

sin i? r A sin w 

so that on 2-forms * 2 = — 1, as should be in accordance with (A. 4). More details about 
the Hodge star operator can be found in |19| . 

At last, it should be noted that all the above is easily over linearity continued on 
the matrix- valued differential forms (see, e. g., Ref. [2T|), i. e., on the arbitrary linear 
combinations of forms di^..^ dx l1 A dx 12 A...Adx lp , where coefficients £i£ 1 i 2 ...i J1 belong 
to some space of matrices V, for example, a SU(3)-Lie algebra. But now the Cartan's 
wedge (exterior) product A should be understood as product of matrices with elements 
consisting of usual (scalar) differential forms. In the SU(3)-case, if T a are matrices of 
generators of the SU(3)-Lie algebra in 3-dimensional representation, we continue the 
above scalar product G on the SU(3)-Lie algebra valued 1-forms A = A^Tadx^ 1 and 
B = B h v T h dx v by the relation 

G(A, B) = g^A«B b u Tr(T a T b ) , (A. 10) 

where Tr signifies the trace of a matrix, and, on linearity with the help of (A.l), G can 
be continued over any SU(3)-Lie algebra valued forms. Such a matrix-valued gener- 
alization of *-operator is extremely useful when exploring solutions of the Yang-Mills 
(and Maxwell) equations in Minkowski spacetime and was systematically employed in 
Refs. PH3]. 
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Operations of vector analysis in terms of * and d 

The basic property of the exterior differential d is the same form in the arbitrary 
curvilinear coordinates x^ on Q. Namely, d = d Xi dx-i with d Xi = d/dxi- For example, 

i 

in Cartesian and spherical coordinates we have, respectively, d — d x dx + d y dy + d z dz 
or d — d r dr + dtfd$ + dpdip. 

Passing on now to the vector analysis on Q, we should note that it is usually formu- 
lated in Cartesian coordinates, where all main operations ( divergence, curl operator 
and so on) look in the simplest form which makes difficulties when transfering to the 
arbitrary curvilinear coordinates. We can, however, simplify the situation if noting 
that there is one-to-one correspondence between any vector field a = a x \ + a-yj + a z \i 
with a x , y ,z = a x ,y,z{x,y, z) and 1-form a x dx + a y dy + a z dz so that the latter will 
be denoted by the same notation a in what follows. When describing vector fields by 
1-forms, however, we at once gain a number of advantages. Indeed, according to the 
standard rules |21| it is easy to get an expression of a in the arbitrary curvilinear coor- 
dinates x l , i = 1,2,3, if knowing x,y,z as the functions x l . This is done by replacing 
dx = (dx I dx l )dx l , dy — (dy/dx l )dx l , dz — (dz / dx l )dx l , x,y,z — x,y,z(x l ) in the 
expression of a in Cartesian coordinates and we at one blow obtain the components 
of a in the given curvilinear coordinates - those are just the coefficients at the cor- 
responding dx % . Further, the basic products of vectors, scalar and vectorial ones, are 
now simply written through the conforming differential forms. Namely, scalar product 
(a, b) = ab = g^a^bu, vectorial product axb = *(a A b), where metric coefficients 
and the Hodge star operator * are defined by the expression ds 2 in the given 
curvilinear coordinates [see, e.g., (A.6)-(A.8)]. 

At last, we have the standard de Rham complex [181119] 

F(f2) = A° -A A 1 A A 2 A 3 -A , (All) 

which enables us to write down the main operations of vector analysis as 

Qf . 

grad f = df = -—^dx % , rot a = *(da) , div a = *(d * a) (A. 12) 

with arbitrary function / 6 F{Q), so that, in virtue of the famous property d 2 = 
for operator d, we automatically obtain the identities rot grad/ = 0, div rot a = 0, 
provided that the first and second de Rham cohomology groups of SI are equal to 
zero: H O = H 2 S1 = 0. It should be emphasized that relations (A. 12) hold true for 
any curvilinear coordinates as soon as the expression of metric ds 2 is fixed in those 
coordinates and, accordingly, the Hodge star operator action is defined on any p-form. 
After computing with using the given curvilinear coordinates we can always return to 
Cartesian ones by replacing dx 1 = [dx 1 / dx)dx + (dx 1 / dy)dy + (dx 1 / dz)dz, dx — > i, 
dy — > j, dz — > k. 

Many relations of vector analysis obtained in Cartesian coordinates can easily be 
generalized to the arbitrary curvilinear coordinates within the formulation under con- 
sideration. For example, the identities 

div(a x b) = (rot a)b - a(rot b) , (A. 13) 



rot (a x b) = (bV)a - (aV)b + div(b)a - div(a)b (A. 14) 
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acquire the form 

div(a x b) = *d(a A b) , (A. 15) 

rot (a x b) = *d * (a A b) (A. 16) 

holding true for any curvilinear coordinates. On the other hand, if trying to use (A. 13)- 
(A.14) within framework of the standard formulation, e.g., in spherical coordinates, 
then this will lead to the perfectly bulky expressions. At the same time, the right-hand 
sides of (A.15)-(A.16) are easily computed for concrete a and b in any curvilinear 
coordinates along the lines above. 



The confining magnetic field 

To illustrate some of the above let us compute the strength H of the confining magnetic 
field of (8)-(9) in Cartesian coordinates. By definition we have H = rot A and in 
spherical coordinates A = (br + B)dip. Then in accordance with (A. 12) and (A. 7) 
H = *(dA) = *d[(br + B)dp>] = *(bdr A dtp) = -g^dtf = H®dd which entails the 



r sin 



t) ■ 



module H = s f^ 7 H^H v = Jg m H$ 

On the other hand, $ = arccos - = arccos . z = wherefrom 

y/x 2 +y 2 +z 2 

, . d$ . d-d , d-d , xz , yz , J x 2 + „ , 

air = ^dx + — dy + — dz = dx H dy = dz 

dx dy dz r 2 ^J x 2 + y 2 r 2^/ x 2 +y 2 r 2 - 
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At last, sin i? = — — and we get 

TT b , . b I xz yz , , 

H = — : — -av — — -= ^dx H = Tjdy — dz 

smv r \x z + y z x A + y A 

or, using the above isomorphism dx <=>■ i, dy <=>■ j, dz <=>- k, 

b I xz yz . , \ , . _ 

H = — i+ " 2 J-k , (A.17 

r \x z + y z x z + y z I 

so the module H = J Hi + H$ + Hi = . |b| = -IJL. 

V x y z ^/ x 2_^_y2 r sin # 



Appendix B 

The facts adduced here have been obained in Refs. [2,3 and we concisely give them 
only for completeness of discussion in Section 2. 

For the case of U(l)-group the Yang-Mills equations (1) turn into the second pair 
of Maxwell equations 

d*F = (B.l) 

with F — dA, A = A^dx^ . The most general ansatz for a spherically symmetric 
solution is A = A t {r)dt + A r (r)dr + A^{r)dd + A v (r)dip. 

For the latter ansatz we have F — dA — —d r AtdtAdr + d r A^drAd'& + drA l pdrAdifi 
for an arbitrary A r (r). Then, according to (A. 9), we obtain 

9 1 
*F = (r sintf)<9 r A t dtf A dip + smfidrA^dt A dip ; — -d r A v dt A d$ (B.2) 
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which entails 

d* F = sin $d r (r 2 d r A t ) dr A d$ A dip - sin ■dd 2 A$dt Adr A dtp- 

1 9 

cos -ddrA-adt A d-d A dip H -dt-A^dt A dr A d$ , (B.3) 

sin if 

wherefrom one can conclude that 

a r (r 2 9 r A f ) = 0, d 2 = , (B.4) 

d?A# = d r A# = . (B.5) 

and we draw the conclusion that A# — C\ with some constant C\ . But then the Lorentz 
condition (2) for the given ansatz gives rise to 

smdd r {r 2 A r ) + d$(sm'&A#) = smfid r {r 2 A r ) + d^d sun?) = 0, 

or 

d r (r 2 A r ) + d cot & = 0, (B.6) 

which yields A r = —C\ cot i? jr + C2/r 2 with a constant G%. But the confining solutions 
should be spherically symmetric and contain only the components which are Coulomb- 
like or linear in r, so one should put C\ = C2 = 0. Consequently, the ansatz A — 
At(r)dt + A v (r)dip is the most general spherically symmetric one and then equations 
(B.4) give 

A t = - + A,A V = br + B (B.7) 
r 

with some constants a, b, A, B parametrizing solutions which proves the uniqueness 
theorem of Section 2 for U(l)-group. Minor modification of the above considerations 
allows us to spread the proof to the Yang-Mills equations (1) (for more details see Refs. 

EE])- 
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